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We investigate how hybridization (single-quasiparticle scattering) between two superconducting 
bands induces odd-frequency superconductivity in a multiband superconductor. An explicit deriva¬ 
tion of the odd-frequency pairing correlation and its full frequency dependence is given. We find 
that the density of states is modified, from the sum of the two BCS spectra, at higher energies by 
additional hybridization gaps with strong coherence peaks when odd-frequency pairing is present. 


Many materials have multiple bands close to the Fermi 
level. It is then not surprising that there also exist many 
superconductors with more than one superconducting 
band. Well-known multiband superconductors are MgB 2 
[l|-Q, which hosts two distinct superconducting gaps, and 
the iron-based superconductors 0-@|, where the order 
parameter changes sign between different bands [7, gj. 
Multiband superconductivity has also been suggested to 
be important in simple metals ‘IfTiT. heavy fermion com¬ 
pounds jll-15], different carbides (Ty, [17|, and Chevrel 
phases |l8[, as well as for engineering time-reversal in¬ 
variant topological superconducting states 3, 20| . 

Multiple superconducting bands allow for unusual cou¬ 
pling effects. Early, and much current, theoretical focus 
has been oriented towards studying the effects of Joseph- 
son coupling between different superconducting bands, 
i.e. the exchange of whole Cooper pairs between bands 
HHH. lowever, conceptually much simpler is single¬ 


quasiparticle scattering, or tunneling, between the su¬ 
perconducting bands. The origin of interband quasipar¬ 
ticle scattering can be impurity scattering, but a com¬ 
mon intrinsic source is superconductivity associated with 
specific orbitals, which subsequently hybridize to form 
the low-energy bands around the Fermi surface. Tun¬ 
neling spectroscopy has found significant effects of inter¬ 
band quasiparticle scattering in silicon clathrate BasSLifi 
|22, .28!], iron-based superconductors [2flj |. and MgB 2 [30j . 
Recent ARPES data on MgB 2 have also shown that in¬ 
terband scattering due to disorder can be important [311. 

Theoretically, interband single-quasiparticle scattering 
results in both band hybridization and interband pairing, 
where the two electrons forming a Cooper pair belong 
to different bands. Straightforward effects of band hy¬ 
bridization were studied already quite early on 32l-34|. 
More recently, band hybridization has also been proposed 
to influence the nodal structure of the superconducting 
gap in iron-based superconductors flB-liist. On the other 
hand, consequences of the induced interband pairing have 
been much less discussed. Pure interband pairing in cold 
atoms and quantum chromodynamics systems has been 
proposed to give a “breached” reg ime containing both a 
superfluid and a normal liquid [39j,|40(, however, in band- 


hybridized superconductors the interband pairing is also 
accompanied by (usually large) conventional intraband 
pairings. Still, it was recently pointed out, using symme¬ 
try arguments and simple mean-field BCS calculations, 
that odd-frequency, odd-interband pairing can be ubiq¬ 
uitous in multiband superconductors 41]. 

The fermionic nature of the superconducting wave 
function usually renders a division into spin-singlet even- 
parity (i.e. s-, d- wave) or spin-triplet odd-parity ( p- 
wave) superconductors, but it can also be even or odd 
under time, or equivalently frequency 4^-44|. Exam¬ 
ples of this are odd-frequency spin-triplet s-wave super¬ 
conductivity giving rise to long-scale proximity effects 
in superconductor-ferromagnet systems 4^-43] or odd- 
frequency spin-singlet p -wave pairing in non-magnetic 


junctions [48H50|. In multiband superconductors the 


band index offers the additional possibility of spin-singlet 
s-wave pairing which is odd in both band and frequency, 
without translation or spin rotation symmetry breaking. 


In this work we derive the exact Green’s functions 
for a generic multiband superconductor with single¬ 
quasiparticle scattering. We thus obtain the full fre¬ 
quency dependence of the interband pairing and find that 
the odd-interband pairing has odd frequency dependence, 
while the even-interband pairing is even in frequency. By 
studying the density of states (DOS) we also discover that 
finite band hybridization gives rise to an extra gap lo¬ 
cated beyond the original gap edges. This hybridization- 
induced gap is not fully depleted, but has very pro¬ 
nounced BCS-like coherence peaks. Moreover, the hy¬ 
bridization gap disappears whenever the odd-frequency 
interband pairing is zero and is thus a directly measur¬ 
able signal of odd-frequency superconductivity. 


To model a generic multiband superconductor we con¬ 
sider two bands, band 1 and 2, with dispersions tki and 
£fc 2 , where e*, = £_&. For simplicity, we assume that 
the two bands independently develop conventional spin- 
singlet s-wave superconducting order parameters Ai and 
A 2 , respectively, but one of them can also be zero. Fi¬ 
nally, we add a small single-quasiparticle hybridization, 
or scattering, term proportional to T between the two 
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bands, resulting in the Hamiltonian: 


r < (w 2 + B 2 ) 4 / 4 (w 2 + -Bl) 1 / 4 , which translates into: 


H = E tklQ'ka-Q'kcr ^k2^\ a ^ka + E r ( fc ) a ^ 6fcCT+H - c - 


kcr 


ka 


+ E A i(*)4t°-fc4. + A ^ k ) b U b -H + Hc -> 


a) 


where a ( b) is the annihilation operator in band 1 (2). 
Alternatively, the band hybridization can be interpreted 
as a coupling process in real space, with the a- and b- 
electrons living to the left and the right of a junction or 
in different layers 32, 5l|. In this picture the Josephson 
coupling would instead be a two-particle tunneling term 
(not included here). 

We start by calculating the spin-singlet s-wave inter¬ 
band anomalous Green’s functions El 2 and F 21 , which 
express the pairing correlations of two electrons belong¬ 
ing to different bands. Assuming T to be a small pa¬ 
rameter, we can use standard perturbation theory |52|. 
The first order contributions are then represented by the 
schematic diagrams in the inset of Fig. [TJa) and give: 
F$ = F 1 TG 2 — ^iTF 2 , where t? is the hole propaga¬ 
tor, i.e., G = —G(—fc, —w). The minus sign before the 
second term is due to scattering of hole propagators (left¬ 
going arrows). The normal and anomalous propagators 
without the hybridization are as usual [52| : 


Gj Fj 

A’ 


M 2 - B 2 


iu + e k j 

A* 


• ), ( 2 ) 

IU - £ k jJ 


where B 2 = E kj = e 2 k - + |A,,-| 2 and u = u n = n{2n + 
1 )ksT are the fermionic Matsubara frequencies. Using 
these expressions we arrive at F$ = r[zw(Ai — A 2 ) + 
Aiefc 2 + A 2 efei]/[(w 2 -|-B 2 )(w 2 -(-B|)]. The next to leading 
order terms for U 12 are cubic in T. Further organizing 
the perturbation expansion of the interband pairing of a 
given order n in a systematic way, we find several recur¬ 
sion relationships 61]. These can be compactly written 
in a matrix form as: 





( - r 12 
p(n — 2 ) 

W^12 


( 3 ) 


where g = F 2 /[(u; 2 + E? 1 ){u 2 + B? 2 )], e = {iu — £ 1 ){iu — 
e 2 ) - AiA 2 and / = -iw( A{ - A 2 ) + Aje 2 + A^. For 
the starting point of the recursion we use the first or¬ 
der anomalous interband Green’s function U 12 and the 
corresponding normal propagator: G ^ = F[AiA 2 — 
{iu — £\){iu — e 2 )]/[(w 2 + E 2 ){u 2 + E 2 )}. From Eq. ([3]) 
we recognize that the Green’s functions to infinite or¬ 
der can be written as a geometric series, where the quo¬ 
tient is a two-by-two matrix. The formal criterion for 
a matrix geometric series to be convergent is that the 
norm (i.e. the largest singular value) of the coefficient 
matrix is < 1. From this we arrive at the condition 


r < V|Ai||A 2 |, (4) 

meaning that the series is always convergent for suffi¬ 
ciently small r, provided that Ai and A 2 are both finite. 
This allows us to sum the infinite series and we arrive 
at F \2 = F[zaz(Ai — A 2 ) + A 2 Ci H- A^e 2 ]/B, where D == 
(w 2 +B 2 )(w 2 +B|)—F 2 [2(e 1 e 2 —w 2 ) —A 2 A!—A|A 2 ]+r 4 . 
The expression for B 2 i is obtained by exchanging the 
band indices and we can also form the odd and even 
combinations of B 12 and B 2 i with respect to the band 
index: 

E° 2 dd (k, ku) = Fl2 2 F21 = iuT{A, - A 2 )ID (5) 

Bf 2 ven (k, iu) = Fl2 F F21 = F(A iefc2 + A 2 £ kl )/D. ( 6 ) 

The odd-band combination is directly seen to be odd in 
frequency, whereas the even-band combination has a con¬ 
ventional even-frequency dependence. This is fully con¬ 
sistent with Fermi-Dirac statistics for spin-singlet s-wave 
superconducting pairing. Furthermore, we see that in¬ 
terband pairing always requires a finite band hybridiza¬ 
tion and that the odd-interband pairing also requires 
Ai 7 ^ A 2 . Equations ((SJ)-® can be Fourier transformed 
to real space and then evaluated numerically, with the 
result shown in Fig. [T] By definition, the odd-frequency 
pairing amplitude must be zero at u = 0. Close to u = 0, 
we find that the leading term in B odd is linear in u. How¬ 
ever, the slope can initially be large and then abruptly 
change sign to asymptotically go to zero for large u, re¬ 
sulting in an approximate l/w-dependence, which has 
also been found for some odd-frequency states [H 3 , 0 ] ■ 

By forming an analogue of Eq. © for the intraband 
anomalous Green’s function in band 1 we find F\ = 
{Ai[(iu ;) 2 — B|] — F 2 A 2 }/B and the corresponding nor¬ 
mal Green’s function: Gi = {{iu + £i)[(iuj) 2 — Bf] — 
T 2 {iu — e 2 )}/B. The expressions for B 2 and G 2 are ob¬ 
tained by mutually exchanging band indices. The ex¬ 
pressions for the normal Green’s functions allow us to 
compute the DOS using the standard formula: 

N(E) = — —ImTr G(B + id), d-»0 + (7) 

7T 

Here the trace involves a sum over band and spin indices 
and an integral over fc-space. We use similar expressions, 
but unsummed over the band index, to define the partial 
DOS N\ and N 2 . 

The total and partial DOS offer a direct connection 
to experimental measurements on multiband supercon¬ 
ductors. In Figs. [2] - [4] we show numerically obtained 
results for the DOS, explicitly exploring the effect of in¬ 
terband pairing. To most clearly illustrate the effect of 
interband pairing we use two generic parabolic bands: 
£j = h 2 k 2 /2irij — /ij, with effective masses mi = 20 m e , 
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Figure 1: Odd- (a) and even- (b) interband pairing amplitudes 
in meV per nm 3 when A 2 = 2.5 meV, Ai = 2.5, 2.8, 4.5 meV 
(blue, orange, green), and T = 3 meV, with the band structure 
specified in the main text. Inset: first order hybridization 
contributions to the interband pairing; F 1 YG 2 (top), — G 1 TF 2 
(bottom). Solid (dashed) line represents the propagator in 
band 1 (2) and x the hybridization. 


TO 2 = 22 m e and distances from the bottom of the bands 
to the Fermi level /.ti = 100 meV and fi 2 = 105 meV. For 
this and other band structures we studied, the interband 
effects are most clearly visible when T is comparable to 
|ei — £ 2 ! for fixed k ~ kp. The contributions to the DOS 
are obtained by numerical integration in the range includ¬ 
ing both Fermi surfaces and we use a smearing parameter 
S = 0.01 meV. We have also independently verified the 
DOS results by exact diagonalization of the Hamiltonian 
in Eq. ©■ In fact, we find a perfect agreement even 
far beyond the theoretical condition for convergence in 
Eq. (©. 

We here especially showcase that unusual features in 
the DOS are only seen when the conditions T / 0 and 
A 1 ^ A 2 are both satisfied, which are exactly the two 
key criteria for odd-frequency pairing, see Eq. ©• First, 
the DOS without any band hybridization, as shown in 
Fig. ©a), is just a sum of two BCS spectra with energy 
gaps E g i = Ai and E g 2 = A 2 , respectively, as expected. 
However, when we turn on hybridization, see Figs. ©b- 
d), we see extra, very notable, dips in the DOS located 
beyond the original gap edges E g 1 ^. These dips, symmet¬ 
ric around zero energy, clearly resemble superconducting 
gaps with their pronounced BCS-like coherence peaks. 
However, the DOS in the gap regions are not zero, but 
instead equal to the partial DOS at these energies. Still, 
we refer to these features as hybridization-induced gaps. 


(a) 

"l 

> 

IS 


(0 


Ai = 2.5 meV, Ag = 1 meV, r = 0 meV 




(b) 

”| 

5 


(d) 


= 2.5 meV, Ag = 1 meV, r = 0.5 meV 

2000 
1500 
1000 
500 
0 

-20 -15 -10 -5 0 5 10 15 20 

E - E f [meV] 

At = 2.5 meV, Ag = 1 meV, r = 6 meV 




N —— 

N, - 

n 2 

. J 



— j 

- y^j 

* 


rt 






Figure 2: Total ( N ) and partial densities of states (IVi, N 2 ) 
when Ai = 2.5 meV and A 2 = 1 meV for different values of 
T = 0, 0.5, 3, 6 meV (a-d), with the band structure specified 
in the main text. 


The hybridization-induced gaps grow in size and move 
to higher energies for larger band hybridizations, and 
we associate these features with interband superconduct¬ 
ing pairing. We note that beyond these hybridization- 
induced gaps, we see no other distinctive features as¬ 
sociated with interband pairing. Notably, there are no 
zero-energy or subgap states, otherwise often associated 
with odd-frequency pairing 43. EHl. IHol [EBI - IHij . Recent 
works have pointed out that zero-energy states do not 
always accompany odd-frequency pairing Ell, [HIM Hil, 
and odd-frequency, odd-interband pairing provides an¬ 
other example when this is not the case. 


In Fig. [3]we instead fix F = 3 meV, and A 2 = 2.5 meV 
but vary Ai. Distinct hybridization-induced gaps are 
present at energies beyond the original gap edges, in¬ 
dependent on the relative size of the two original gaps. 
The only exception is exactly when Ai = A 2 , then 
the hybridization-induced gaps completely disappear, de¬ 
spite the finite band hybridization, see Fig. 0c). Detun¬ 
ing the value of Ai slightly from that of A 2 results in 
small, but noticeable, dips in the DOS, at the positions 
where the full gaps develop for increasing differences be¬ 
tween Aiand A 2 . We thus find that the hybridization- 
induced gaps are only present in the DOS when both 
r / 0 and Ai ^ A 2 . These are exactly the two key 
criteria for odd-frequency pairing, as seen in Eq. ©■ 
In fact, only the odd-frequency interband pairing dis¬ 
appears at Ai = A 2 , the even-frequency part is in gen¬ 
eral non-zero as soon as T / 0. Specifically, the even- 
and odd-frequency interband pairing amplitudes corre¬ 
sponding to the parameters in Fig. ©c)-(e) are plotted 
in Fig. [T] From there it is clear that the even-frequency 
interband pairing is large and not changing significantly 
around ui = 0 , while the odd-frequency part changes 
from identically zero in Fig. ©c) to a notable non-zero 
derivative at w = 0 for Figs. ©d )-( e). We can thus con- 
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A 1 = 0.5 meV, = 2.5 meV, r = 3 meV 






A 1 = 2 meV, Ag = 2.5 meV, r = 3 meV 




A 1 = 7.5 meV, Ag = 2.5 meV, r = 3 meV 



(h) A, =15 meV, A 2 = 2.5meV, r = 3 meV 



Figure 3: Total (N) and partial densities of states (Ah, N%) 
when A 2 = 2.5 meV and T = 3 meV for different values of 
Ai = 0.5, 2, 2.5, 2.8, 4.5, 7.5,10, 15 meV (a-h), with the band 
structure specified in the main text. 


elude that odd-frequency interband pairing is necessary 
for producing the hybridization-induced gaps. Detect¬ 
ing gaps beyond the original two superconducting gaps 
in multiband superconductors is therefore a clear sign of 
the presence of odd-frequency pairing. Intriguingly, ad¬ 
ditional gap features have already been reported in the 
multiband superconductor BasSLi6 [271 . 

Finally, we also show that only one band has to be 
natively superconducting for the hybridization-induced 
gaps to be present. In Fig. 0] we display how the 
hybridization-induced gap grows with increasing T when 
A 2 = 0. Finite single-particle hybridization T results in 
a proximity-induced gap also in the second band, which 
is manifested as a gap around zero energy, although al¬ 
ways smaller than E g 1 = Ai. For all parameter choices 
in Fig. [I] the spectrum always has a gap at zero energy, 
i.e. there are no zero-energy states. This is not clearly 
seen in Fig. |4ja) due to the finite smearing parameter, 
but it is clearly visible in exact diagonalization results 
which do not suffer from the same problem. In addition, 
a finite T results in finite odd-frequency interband pair¬ 
ing, and we consequently also see hybridization-induced 


(a) 


A| = 3 meV, Ag = 0 meV, r = 1 meV 



(c) 


A| = 3 meV, Ag = 0 meV, r = 3 meV 



( b ) A, = 3 meV, Ag = 0 meV, r = 2 meV 




Figure 4: Total ( N ) and partial densities of states (Ah, Ah) 
when Ai = 3 meV and A 2 = 0, i.e. only one band with native 
superconductivity, for different values of T = 1, 2, 3, 5 meV (a- 
d) with the band structure is specified in the main text. 


gaps beyond the Ai gap, which also grow with F. In 
fact, these gaps are much more pronounced than the 
proximity-induced gap in band 2 at zero energy. 

In summary, we have studied the effect of single¬ 
quasiparticle hybridization or scattering in a two-band 
superconductor. By performing perturbation theory 
to infinite order in the hybridization term, we have 
obtained the exact, fully frequency dependent, ex¬ 
pression for the interband pairing, which can be di¬ 
vided up into odd-frequency odd-interband and even- 
frequency even-interband pairing. The conditions for fi¬ 
nite odd-frequency interband pairing are (a) finite single¬ 
quasiparticle hybridization and (b) a non-zero difference 
between the original superconducting gaps; no applied 
magnetic field, inhomogeneity, or interface is required. 
Furthermore, we have shown that the DOS develops non¬ 
trivial gaps features with distinct coherence peaks be¬ 
yond the original gap edges only if the conditions for 
odd-frequency pairing are satisfied, otherwise the spec¬ 
trum just a sum of two BCS spectra. Detecting such 
additional gaps thus provides experimental evidence of 
odd-frequency pairing in multiband superconductors. 
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In this Supplementary Material we derive in detail the 
exact result (within perturbation theory) of the inter¬ 
band pairing amplitude in a two-band superconductor 
with hybridized bands. 

The Hamiltonian for a generic multiband superconduc¬ 
tor can be written as (same as Eq. (1) in the main text): 

H = E £kial a CLka + ek2b' ka bkcr + E r ^4A-+H.c. 

key key 

+ E A i(*)4r a -k4. + Mk)b{ t bl H + H.c.. (SI) 

k 

It describes two superconducting bands coupled by a hy¬ 
bridization, or scattering, term T(fc), which we treat as a 
perturbation. We thus start with two copies of the single 
band superconducting Green’s functions (same as Eq. (2) 
in the main text): 


(S2) 


which we use together with the hybridization term to 
build up the full Green’s functions of the system. Here 
G is the normal electron propagator, which we schemat¬ 
ically denote by->-, and G is the hole propaga¬ 
tor denoted by-<-. Moreover, F is the anoma¬ 
lous propagator, denoted by —<->—, and F * is its 

Hermitian conjugate —> -<—. We will furthermore 

use solid lines to denote the propagators in band 1 and 
dashed lines for band 2. The hybridization term is rep¬ 
resented by x, which always connects two propagators 
from different bands with the same direction of the ar¬ 
row (due to momentum conservation). 

Now, we want to calculate the interband pairing am¬ 
plitude F 12 , which corresponds to the sum of all pro¬ 
cesses of the type < C '.A > . The simplest processes, 
here indicated by the superscript (1), include just one 

scattering event x: F^ = —<-X.>. — 

-<- X— <-.>-- = F\ I'Gj — t?iTF 2 . The sign 

of a particular process is simply given by (—1) ; , where 
l is the number of hole scattering events, i.e. scattering 
connecting two left-going arrows. 

Since Fj 2 of any given order has to end with a right¬ 
pointing arrow in band two, there are only two possibili¬ 
ties, either it ends with X .>.or with X <.>--. 

Similarly, it needs to start with a left-going solid arrow, 
which also limits the options. By drawing all possible 
processes and translating them into formulas, we can 
write: 

f[2 - T;\" 1 x <.> - /V" IS X.>.. 

(S3) 


Carrying this procedure one step further, we write 
E n-1) and using f£~ 2) and 

t { r 1] =f 1 ( 2 - 2 ) >—<—Er 2) *—<— 

(54) 

F ( n - 1) =F 1 ( 2 _2) X->-t?£p 2) *-s->—• 

(55) 


Plugging these into Eq. (lS3l) we arrive at: 


rp(n) VT(n—2)r 1 

/ |.j = Cr 12 |x-•-x < • -> x • | 4 - 

f£ 2 2) [x->-x-—•>-— x >-< x <-> j, (S6) 

By a similar procedure we get: 

G12 = G’ • 'j 2) [x-•-x-.<. \ <-X x >. <.] + 

Fy 2 2) [x-•>-x > .<.• \ --- x - ]. (S7) 


Thus, we get a closed set of equations if we con¬ 
sider F \2 and G12 together, which is easiest done in 
a matrix formalism. Note that all arrow diagrams 
can be directly translated to specific formulas using 

Eq. ©, e.g. x-<-x <.> = - l'G^ \')F, = 

( _r) 2 (* w - ei )A 2 /{[M 2 - Ef] [(iw) 2 - El}}. 

Summarizing, we arrive at the matrix recursion rela¬ 
tion in Eq. (3) in the main text: 



= 9 


-r 



(S8) 


where g = T 2 /[(w 2 +£ 2 j)(w 2 + F 2 2 )], e = (iu> — ei)(iw — 
e 2 ) — AiA 2 and / = — iw(AJ — A 2 ) + A*e 2 + A^ei- 
This is a matrix geometric series and is thus easily 
summed giving: 



gf 

1 -ge 



(S9) 


where q = [(1 - ge){ 1 - ge*) + g 2 \f \ 2 ] \ = 

—^-5-X< >.< - -<-X < — = 

T[AiAJ - (iuj - ei)(iu - e 2 )]/[(w 2 + Ef)(uj 2 + El)}. 
Some terms cancel and we get Fj 2 = qF[2 ■ Finally 
F \2 = r[zw(Ai - A 2 ) + A 2 ei + Aie 2 ]/F, where D = 
(w 2 +F 2 )( w 2 +F2 2 )-r 2 [2(e 1 e 2 - W 2 )-A^A 1 -A|A 2 ]+r 4 . 

For the normal propagator Gi, needed for calculating 
the DOS, we create a matrix equation together with E\. 
otherwise the procedure proceeds in the same manner as 
above. 































